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Lattice field theory provides a quantitative tool to study the role of nonpertur- 
bative semiclassical configurations in QCD. This talk briefly reviews our present 
_ understanding of the role of instantons in QCD and describes in detail new devel- 

(_P , opments in the study of morons on the lattice. 

(N 



1 Introduction 
> 

t^^ ■ As highlighted by this conference, despite a quarter century of effort, under- 

standing the nonperturbative origin of confinement in QCD has remained an 
elusive challenge. In this talk, I will describe how combining lattice field the- 

^~^ I ory with familiar semiclassical ideas is providing new insight into the role of 

f^ . instantons and nierons in the QCD vacuum and new tools to understand the 

^D ' physics of confinement. 

One of the great advantages of the path integral formulation of quantum 
mechanics and field theory is the possibility of identifying non-perturbatively 
the stationary configurations that dominate the action and thereby identify 
(Vj " and understand the essential physics of complex systems with many degrees of 

<~i . freedom. Thus, the discovery of instantons in 1975ld gave rise to great excite- 

ment and optimism that they were the key to understanding QCD. Indeed, in 
contrast to other many body systems in which the quanta exchanged between 
interacting fermions can be subsumed into a potential, it appeared that QCD 
S . was fundamentally different, with topological excitations of the gluon field 

dominating the physics and being responsible for a host of novel and important 
effects including the 6 vacuum, the axial anomaly, fermion zero modes, the mass 
of the 77', the chiral condensate, and possibly even confinement. However, until 
lattice QCD became sufficiently sophisticated and sufficient resources could be 
devoted to the study of instanton and meron physics, it had not been possible 
to proceed beyond tha dilute instanton gas approximation and a qualitative 
exploration of meronsH. 

Using lattice field theory, our basic strategy is to reverse the usual analyt- 
ical process of calculating stationary configurations and approximately sum- 
ming the fluctuations around them. Rather, we use Monte Carlo sampling of 
the path integral for QCD on a lattice to identify typical paths contributing to 



the action and then work backwards to identify the smooth classical solutions 
about which these paths are fluctuating. 

2 The Role of Instantons 

Insight into the role of instantpns and their zero modes from lattice QCD has 
been reviewed in detail in refp, and I will only summarize the highlights her£. 
The instanton content of gluon configurations can be extracted by cooling □, 
and the instanton size is consistent with the instanton liquid model and the 
topological susceptibility agrees with the Veneziano-Witten formula B. We 
obtain striking agreement between vacuum correlation functions, ground state 
density-density correlatioa functions, and masses calculated with all gluons 
and with only instantonsH. Zero modes associated with instantons are clearly 
evident in the Dirac spectrum, and account for the p, tt, and 77' contributions 
to vacuum correlation functions U. Finally, we have |Observed directly quark 
localization at instantons in uncooled configurations lI. The physical picture 
that arises corresponds closely to the physical arguments and instanton models 
in which the zero modes associated with instantons produce localized quark 
states,_and quark propagation proceeds primarily by hopping between these 
states El. However, of particular relevance to this conference, it is now clear 
that simple configurations of instantons, such as a dilute gas a or a random 
superposition □, do not confine color charges. It is especially important to 
emphasize that even an ensemble of instantons with a distribution of large 
instantons oc p~^ yields a static potential which approaches a constant at large 
distance when one is sufficiently careful in calculating the contributions of very 
large Wilson loopsil. Although it is possible that instanton configurations that 
differ essentially from a random superposition could lead to confinement, it is 
interesting to discuss here closely related topological gluon excitations known 
as merons. 

3 Merons on the Lattice 

In the space available in these proceedings, I would like to concentrate on new 
advances with James V. Steele u in studying merons on the lattice. Merons 
are solutions to the classical Yang-Mills field equations with one-half unit of 
localized topological charge. Since they are the (3-|-l)-dimensional nonabelian 
extension of the 't Hooft-Polyakov monopole ccpfiguration, which has been 
shown to lead to confinement in 2-t-l dimensions 113, they have long been con- 
sidered a possible mechanism for confinement a. 

A purely analytic study of merons has proven intractable for several rea- 
sons. Unlike instantons, no exact solution exists for more than two meronsliH, 



gauge fields for isolated merons fall off too slowly {A ^ l/x) to superpose 
them, and the field strength is singular. A patched Ansatz configuration that 
removes the singularities a does not satisfy the classical Yang-Mills equations, 
preventing even calculation of gaussian fluctuations around a meron pairE3. 
Because of these analytic impasses, we have constructed smooth, finite action 
configurations corresponding to a meron pair and their associated fermion zero 
modes on the lattice. For simplicity, we restrict our discussion to SU(2) color 
without loss of generality. 

Two known solutions to the classical Yang-Mills equations in four Eu- 
clidean dimensions are instantons and merons. Both have topological charge, 
can be interpreted as tunneling solutions, and can be written in the general 
form (for the covariant derivative 13^ ^ d^ — iA'^a"- /2) 

Alix) = ^^^/(x^) , (1) 



"^"^ for an instanton and /(x^) — - for a meron. 



with J{x ) ~ X /{x'- + p I iwi c^ii iixo>,c.xi>,^^xi coix^ J y^ ) — 2 

Conformal symmetry of the classical Yang-Mills action, in particular under 
inversion x^ — > x^/x"^ , shows that in addition to a meron at the origin, there is 
a second meron at infinity. These two merons can be mapped to arbitrary po- 
sitions, which we define to be the origin and d^. After a gajjge transformation, 
the gauge field for the two merons takes the simple formllj 

" [x-df 



Ki^)=Va, 



(x - df 



(2) 



Similar to instantons 13, a meron pair can be expressed in singular gauge by 
performing a large gauge transformation about the mid-point of the pair, re- 



sulting in a gauge field that falls off faster at large distances {A 
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caxeful treatment of the singularities shows that the topological charge density 



Q{x) 



Tr 



(f.^F^^'^) ^\5\x) 



1 



+ -r(a;-d) 



(3) 



167r2 V^"' J 2" '"' ' 2 
yielding total topological charge Q — 1, just like the instanton. 

The gauge field Eq. (g) has infinite action density at the singularities x^ — 
{0, d^}. Hence, a finite action Ansatz has been suggestedB 
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^ul-^) — Vafj.iyX < 



X^ 
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x"^ < r , 
r < \fx? < R , 
R< Vx^ . 



(4) 




Figure 1: Meron pair separated by d ■- 
i topological charge. 



/Rr regulated with instanton caps, each containing 



Here, the singular meron fields for v x^ < r and vx^ > R are replaced by 
instanton caps, each containing topological charge ^ to agree with Eq. (||). 
The radii r and R are arbitrary and we will refer to the three separate regions 
defined in Eq. (0) as regions I, II, and III, respectively. The action for this 



configuration is S = ^^ + ^^ In — , which shows the divergence in the r — > 
or i? ^ cx) limit. There is no angular dependence in this patching, and so the 
conformal symmetry of the meron pair is retained. Although this patching of 
instanton caps is continuous, the derivatives are not, and so the equations of 
motion are violated at the boundaries of the regions in Eq. (H). 

Applying the same transformations used to attain Eq. (pp^to the instanton 
cap solution, regions I and III become four-dimensional spheres each containing 
half an instanton 0. The geometry of the instanton caps is shown in Fig. |l| for 
the symmetric choice of displacement d = yRJr along the z-direction. Note 
that the action for the complicated geometry of Fig. |l| with d = ^/Rt can be 
rewritten as 

5 = ^ofl + ^ln-^ , (5) 



with 5*0 = 87r^/(7^. We will concentrate on this case below, and generalization 
to a different d is straightforward. 

The original positions of the two merons Xfj, — {0, dp,} are not the centers 
of the spheres, nor are they the positions of maximum action density, which 
occurs within the spheres with ^max = (48/5^)(i? + r)^/d® at 



(Ml 



^2 '^M . 



(Miii)„ = 



i?^ 



i?2 



d^'^- 



(6) 



However, in the limit r —^ and i? — > oo holding d fixed, the spheres will shrink 
around the original points reducing to the bare meron pair in Eq. (y). In the 
opposite limit R —i- r, the radii of the spheres increase to infinity, leaving an 
instanton of size p = d. 

An instanton can therefore be interpreted as consisting of a meron pair. 
This complements the fact that instantons and antiinstantons have dipole in- 
teractions between each other. If there exists a regime in QCD where meron 
entropy contributes more to the free energy than the logarithmic potential 



between pairs, like the KosteiJitz-Thouless phase transition, instantons could 
break apart into meron pairs 0. The intrinsic size of the instanton caps would 
be determined by the original instanton scale p. 

The Atiyah-Singer index theorem states that a fermion in the presence of 
a gauge field with topological charge Q, described by the equation 

0^ = A^, with ij= (^A , (7) 

has nn right-handed and ul left-handed zero modes (defined by A = 0) such 
that Q =p-pL — nji. This statement can often be strengthened by a vanishing 
theoremlla. Applying p twice to ip decouples the right- and left-hand compo- 
nents, and focusing on the equation for ipn, gives 

D^ + lva^.'J''F^^A Vfl = , (8) 

where fja^u denotes the 't Hooft symbol, cr'^ acts on the spin indices of rpj^,, and 
F^" acts on the color indices. For a self-dual gauge field (like an instanton), the 
second term in Eq. (jg) is zero; and since D^ is a negative definite operator, there 
are no normalizable right-handed zero modes, implying Q — ul. Although the 
meron pair is not self-dual, the second term can be shown to be negative 
definite as well, leading to the same conclusion. 

In general, a gauge field in Lorentz gauge with Q — 1 can be written in 
the form 

A^ix) = -7^a^..d,lnUix) . (9) 

This has a fermion zero mode given by 

i^=(^?\ : C=AAn3/2£-, (10) 



with normalization J\f and e — i<J2 coupling the spin index a to the color index 
a (both of which can be either 1 or 2) in a singlet configuration. The gauge 
field for the meron pair with instanton caps corresponds to 



26^2 



x^+S,f[x-dY 



, for regions J = I, III, 



n(x) ^ { ^^ (11) 

, region II, 

^/x^ix — (i)2 
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Figure 2; The action density S{x) for the regulated meron pair, normalized by the instan- 
ton action So = 8n-^ /g'^, sliced through the center of the configuration in the direction of 
separation. Shown are the analytic (solid), initial lattice (open circles) and cooled lattice 
(filled circles) configurations. The same is shown for the topological charge density Q{x) and 
fermion zero mode density ip'^^l;{x). 




with ^i = rjd and ^m 
then Eq. ^ with 



Rj d. The normahzed solution to the zero mode is 



M-^ = 2TT(f 



2- W- 



Note that the uaxegulated meron pair (r 
zero mode itselfll3 



1/2 



0, R 



del 



2\/27r(a;2(x-d)2) 



2^3/4 



(12) 

cxd) has a normalizable 

(13) 



The gauge- invariant zero mode density ^^^(x) has a bridge between two 
merons that falls off like x^^ in contrast to the x~^ fall-off in all other di- 
rections. This behavior can be used to identify merons when analyzing their 
zero modes on the lattice, similar to what was done for instantons in Ref. Q. 

As mentioned above, the patching of instanton caps to obtain an explicit 
analytic solution has unavoidable and unphysical discontinuities in the action 
density. Therefore, in order to study this solution further, we put the gauge 
field on a space-time lattice of spacing a in a box of size Lq x Li x L2 x L3. 






Figure 3: The topological charge density Q{x) on the left and fermion zero mode density 
i/'t'!/>(x) on the right for a cooled meron pair with r = 9, separated by distance d = 21, 18, 
and 14 respectively, in the {z, i)-plane. 



The gauge-field degrees of freedom are replaced by the usual parallel trans- 



port, U^{x) — Pexp 



The exponentiated integral can 



be performed analytically within the instanton caps, producing arctangents. 
Outside of the caps, the integral is evaluated numerically by dividing each link 
into as many sub-links as necessary to reduce the 0{a^) path-ordering errors 
below machine precision. 

We then calculate the action densit3pS'(2:) for a meron pair with instanton 
caps using the improved action of Ref. E3 and the topological charge density 
Q{x) using products of clovers. These (open circles) are compared with the 
patched Ansatz results (solid line) in Fig. g for r — 9 and d = 21 (in units of 
the lattice spacing) on a 32'^ x 40 lattice. We use the Arnoldi method to solve 
for the zero mode of this gauge configuration and hence the density -0^-0(0;), 
which is also compared with the analytic result in the same figure, showing 
excellent agreement. 



Two important features of the patched Ansatz also evident with the lattice 
representation are the discontinuities in the action density at the boundary of 
the instanton caps and vanishing of the topological charge density outside of 
the caps as given by Eq. (||) . The discontinuity in the action density is unphys- 
ical and can be smoothed out by using a relaxation algorithm to iteratively 
minimize the lattice actional. On a sweep through the lattice, referred to as a 
cooling step, each link is chosen to locally minimize the action density. Since a 
single instanton is already a minimum of the action, this algorithm would leave 
the instanton unchanged (for a suitably improved lattice action). The result 
for the regulated meron pair after ten cooling steps is represented in Fig. g by 
the filled circles, showing the discontinuities in the action density have already 
been smoothed out. In Fig. ^, we also plot side- by-side the cooled Q{x) and 
ip^'ipix) in the (z, i)-plane for three different meron pair separations. As the 
separation vanishes, the zero mode goes over into the well-known instanton 
zero mode result. 

Like instanton-antiinstanton pairs, however, a meron pair is not a strict 
minimum of the action, since it has a weak attractive interaction Eq. (m and 
under repeated relaxation will coalesce to form an instanton. Nevertheless, 
just as it is important to sum all the quasi-stationary instanton-antiinstanton 
configurations to obtain essential nonperturbative physics, meron pairs may be 
expected to play an analogous role. A precise framework for including quasi- 
stationary meron pairs is to introduce a constraint Q[A] on a suitable collective 
variable q (which we chose to be the quadrupole moment 3z^ — x^ — y'^ — t^ of 
the topological charge) as follows 

Z = [da exp{-S[A]} = fdq Zq , (14) 

Z, = [da exp [-S[A] ~ X {Q[A] - qf} . (15) 

The meron pair is then a true minimum of the effective action with constraint, 
allowing for a semiclassical treatment of Zq. Afterwards, q is integrated to 
obtain the full partition function Z. 

The criterion for a good collective variable q of the system is that the 
gradient in the direction of q is small compared to the curvature associated 
with all the quantum fluctuations. In this case, we have an adiabatic limit 
in which relaxation of the unconstrained meron pair slowly evolves through a 
sequence of quasi-stationary solutions, each of which is close to a corresponding 
stationary constrained solution. Detailed comparison of the quasi-stationary 
and constrained solutions shows this adiabatic limit is well satisfied. Therefore, 
we will present here the action of a meron pair as it freely coalesces into an 
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Figure 4: Action of a regulated meron pair as 
tions of dinit = 12, 14, 18, 22 and rinit = 4 on 
trajectories for each cooling step. 

instanton. To compare with the patched 
d and radii r and i? of the lattice configur 
the separation of the two maxima in the 
cubic sphnes) and their values (which 
denoted by S'max — 48/g^w;^), and then 



a function of d/r for initial lattice configura- 
a 32"^ X 40 lattice. The circles show cooling 

Ansatz, we need the meron separation 
ation. We find these by first measuring 
action density A = \M\ — Mni| (using 
are the same in the symmetric case, 
using Eq. ^ to give 



In Fig. 0, we plot the total action of a meron pair as a function of d/r (solid 
line), which for the analytic case is given by Eq. (|g). Also shown are cooling 
trajectories for four lattice configurations with different initial meron pair sep- 
arations. Each case starts with a patched Ansatz of a given separation and the 
lattice action agrees with that of the analytic Ansatz. 

The first few cooling steps primarily decrease the action without changing 
the collective variable (which is now effectively d/r), as observed above in 
Fig. ||. Further cooling gradually decreases the collective variable, tracing out 
a new logarithmic curve for the total action (dashed line in Fig. 4), which 
is about 0.25 lower than the analytic case. This curve represents the total 
action of the smooth adiabatic or constrained lattice meron pair solutions. 
The essential property of merons that could allow them to dominate the path 
integral and confine color charge is the logarithmic interaction Eqs. (pf) which is 
weak enough to be dominated by the meron entropyH. Hence, the key physical 



result from Fig. is the fact that smooth adiabatic or constrained lattice meron 
pair solutions clearly exhibit this logarithmic behavior. 

In summary, we have found stationary meron pair solutions, shown that 
they have a logarithmic interaction, and shown that they have characteristic 
fermion zero modes localized about the individual merons. Hence the presence 
and role of merons in numerical evaluations of the QCD path integral should 
be investigated, and the study of fermionic zero modes and cooling of gauge 
configurations provide possible tools to do so. 
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